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Abstract 

Traffic jam in an optimal velocity model with the backward reference function 
is analyzed. An analytic scaling solution is presented near the critical point 
of the phase separation. The validity of the solution has been confirmed from 
the comparison with the simulation of the model. 
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Recently, the importance of cooperative behavior in dissipative systems consisting of 
discrete elements has been recognized among physists. As a result, granular materials have 
been studied extensively from physical point of views. Similarly, to know the properties 
of traffic jams in daily life is also an attractive subject not only for engineers but also 
physists [01 . There are some similarities between two phenomena in particular in the simplest 
situation where cars and particles are respectively confined in a highway and a long tube. 
Thus, it is interesting to clarify common and universal mathematical structure behind these 
phenomena. 

We propose here a model of the traffic flow 



where Xn and a are the positions of n th car, and the sensitivity, respectively. This model 
contains the psychological effect of drivers. Namely, the driver of x„ takes care of not only 
the distance ahead x„+i — x„ but also the backward distance x„ — The optimal velocity 
function U should be a monotonic increasing function of the distance of Xn+i — Xn and V — 1 
should be a monotonic decreasing function of x„ — Thus, we adopt 



for the later explicit calculation. We put these optimal velocity functions as the product 
form UV in (Q), because the driver of x„ cannot accelerate the car without enough the 
forward distance Xn+i — Xn even when the distance Xn — Xn-i becomes short. This model 




(1) 



U{h) = tanh(/i - 2) + tanh(2); V{h) = 1 + /o(l - tanh(/i - 2)) 



(2) 




(3) 



Our model is also similar to the model of granular flow in a one dimensional tube 




(4) 
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where the exphcit forms of U and the force g are not important in out argument. Although 
real systems contain variety of cars and particles and higher dimensional effects, we believe 
that the most essential parts of both traffic flows and granular pipe flows can be understood 
by pure one dimensional models (|I]) and The reason is as follows: It is known that 
model supplemented by the white noise produces a power law in the frequency spectrum 
of the density correlation function S{q,uj) ~ uj~'^^^, whose exponent 4/3 is very close to 
the experimental value @] and that by the lattice-gas automata simulation [^]. From this 
success the essential effects of randomness such as passing cars and variety of cars seem to 
be represented by the adding white noise to the models (0) and @. 

Komatsu and Sasa reveal that the original OV model can be reduced into the modified 
Korteweg-de Vries (MKdV) equation at the critical point (the averaged car distance h = 
2) for the phase separation. They also show that symmetric kink solitons deformed by 
dissipative corrections describe a bistable phase separation. The exactly solvable models 
in which the essential characteristics of the optimal velocity model are included have been 
proposed ||^. However, as will be shown, the generalized optimal velocity model (|I|) and 
granular model as well as the fluid model of traffic flows by Kerner and Konhauser ^ 
and two fluid models in granular flows ^ are not reduced to MKdV equation but exhibit 
the phase separations between a linearly unstable phase and a stable phase . Thus, there 



is a wider universality class of dissipative particle dynamics which contains (|I|), (H) and fluid 
models 

The aim of this Letter is to obtain an analytic scaled solution of (|l|). To demonstrate 
quantitative validity of our analysis we will compare it with the result of our simulation. 
After the completion of our analysis on (|l]), we will briefly discuss the relation of the result 
and the expected results in and fluid models. 

Let us rewrite (|1]) as 

Tn = a[U{h + rn+i)V{h + r„) - U{h + r„)T/(/i + r„_i) - r„] (5) 
where h is the averaged distance of successive cars and r„ is x„+i — Xn — h . 



Now, let us consider the linear stability of (|^). The linearized equation of (j^) around 
rn{t) = is given by 

rn = a[U'{h)V{h){rn+i - r„) + U{h)V'{h){rn - r„-i) - r„] (6) 

where the prime refers to the differentiation with respect to the argument. With the aid of 
the Fourier transformation rg(t) = Y,n=i exp[— 'ign/i]r„(i(;) with q = 27rm/Nh and the total 
number of cars N we can rewrite (|^) as 

(dt - a4q)){dt - a4q))r,{t) = (7) 

with 



cr±(g) = -2 ± V («/2)' - ^^Dnp, V]{1 - cos{qh)) + ta{UVy sin(g/i), (8) 

where we drop the argument h in U and V. Dh[U,V] = U'{h)V{h) — U{h)V'{h) denotes 
Hirota's derivative. The solution of the initial value problem in (|^ is the linear combination 
of terms in proportion to exp[cr_|_(g)i(:] and exp[cr_(g)t]. The mode in proportion to exp[cr_(g)t] 
can be interpreted as the fast decaying mode, while the term in proportion to exp[cr+(g)i] 
is the slow and more important mode. 

The violation of the linear stability of the uniform solution in (P) is equivalent to 
i?e[o"+(g)] > 0. Assuming qh {qh = is the neutral mode), the instability condition is 
given by 2(UVy'^ cos^( > o-DhlU, V]. Thus, the most unstable mode exists at qh ^ and 
the neutral curve for long wave instability is given by 

The neutral curve in the parameter space (a, h) is shown in Fig.l for /o = 1/(1 + tanh(2)) 
in (0). For later convenience, we write the explicit form of the long wave expansion of cr+ in 
the vicinity of the neutral line 

. .e„,. - cf-^(,l^f - - + 0((,ft)=) (10) 



where cq = (UV)'. Thus, the uniform state becomes unstable due to the negative diffusion 
constant appears for a < an{h). 

The simplest way to describe hydrodynamic mode is the long wave expansion. It is 
easy to derive the KdV equation near the neutral curve from (|I|) as in the case of fluid 
models To describe the phase separations, however, we should choose the critical 

point {a,h) = {ac,hc) from the cross point of {U{h)V{h))" = where the coefficient of 
dxT^ becomes zero and the neutral curve, because KdV equation only has pulses while cubic 
nonlinear terms can produce the interface solution to connect two separated domains. The 
explicit critical point of (0) with /o = 1/(1 + tanh(2)) is given by 

512 

K = 2- tanh~i(l/3) - 1.65343; = /o ^ 1.63866. (11) 

81 

Unfortunately, the reduced equation based on the long wave expansion of our model is an 
ill-posed equation. In fact, the scalings of variables as r„(t) = er{z,T) ,z = e(x + cot) and 



r = e^t with e = ^ ~ a) / Uc leads to 

d^r = aid.r^ - aidr^r + Og^^r^ (12) 

in the lowest order, where oi, 02 and 03 are constants. Its linearized equation around r = do 
is unstable for all scale, because the solution with r — rfo — exp[ikz + A^r] has the growth 
rate Re[Xk] = 2k'^a^dQ which is always positive when 03^0 > 0. 

Of course, this irregularity in the short scale is from the long wave approximation. The 
regularity of the original model (|l]) can be checked easily as follows: Let x = nhhe regarded 
as a continuous variable. From r{x ± h,t) = exp[±hdx]r{x,t) or the Fourier component 
exp[±ig/i] of the translational operator exp[±hdx], the function of the translational operation 
in the shortest scale {qh = vr in the Fourier space) is r{x±h,t) — * —r{x,t). Thus, our model 
in (|^) for r(x,t) in the shortest scale is reduced to 

d^r = a[W{r) - dtr], W{r) = U{h ~ r)V{h + r) - U{h + r)V{h - r). (13) 

Substituting (0) into ([l3| ) it is easy to show W'{r) = — (sech^(/i — r — 2) + sech.^{h + r — 



2)){1 + /o(l + tanh(2))} < 0. Then, the growth rate of the hnearized equation of (|T3|) by 
r — do ~ exp[At + ikx] is given by 



A = ^ ^ (14) 

where Re[X\ < for any d^. Thus, the original model (|l|) is stable for the perturbation in 
the short scale. 

Although it is possible to derive a regularized long wave equation thanks to the Pade 
approximation [11|, the result is more complicated than the original model (|l|). Thus, to 



obtain the asymmetric propagating kink solution, we only eliminate the fast decaying mode 
in (|^) as 

{dt-a+{d,))r{x,t) = (a+-a_)-^iV[r(x,t)], (15) 

where N[r\ represents the nonlinear terms coming from UV . Since ((T+ — cr_)~^ is the inverse 

of the polynomial of the differential operators, it is convenient to use the expansion (cr+ — 
2h 

cr_)^^ ~ a^^[l {UV)'dx + 0{h^)]. Equation (p!5|) is the regularized partial differential 

equation. 

To obtain the scaled propagating kink solution we assume the scaling of the variables by 

e = 



y^(a7--"a)7a^ as 



6c , /6c, X . , 2 



r{x,t) = e^j^jj^^^R{z), z = — {- + cot - e ct) (16) 
where the argument is fixed at h = h^, and c is the positive free parameter which will be 



determined from the perturbation analysis. Substituting ([Tq) into (15) and use the expansion 



N[r]/a = X] E h'^Cmnd^r'^ - h^U'V'd^rdlr + ■■■ (17) 

n=l m=2 

where C21 = \{UV)\C22 = V^]', C23 = ^(t/V^)",C3i = |(f/K)'", C32 = 

j:^Dh[U, V]", C^i = j^{UV)"", and integrate by part of (|lD, we obtain 

"^'^ - R{R^ - 1) + /?^(i?') = e.[^M[R], (18) 



dz"^ dz V Co 



where /? = 3Dh[U, V]' / {2^J co\{UVy"\) . Here we neglect the contribution from the boundary 
and 

,^r^n ( dR\^ dR^ . 1 , dR d^R codR. , 



where l/r^ = v^Z},[f/, V]/co, , P23 = Sv^f/V/ Jco| (f/\/)'"|, P32 = v/3/2D,[f/, \/]7| ([/\/) 



and P41 = V3c^(f^V')"'7(2V2|(t/l^)'"|3). Assuming R{z) = Ro{z) + eRi{z) + ■ ■ we obtain 
an asymmetric kink-antikink solution 

Roiz) = tanh(^±2); e± = ^ ^ ^ ^ (20) 

in the lowest order. The linearized equation of (^) can be reduced to 

CR, = ^M[Ro]; £=^ + l-3i?g + 2/?(i?o^ + ^). (21) 

The solvability condition to determine c is 

/oo 
dz%M[Ro] = 0, (22) 
-00 

where $0 satisfies £^$0 = 0. The explicit form of $0 is given by 

%{z) = (sech[^±2])^/^±. (23) 

Thus, the solvability condition is reduced to 

/ 7-(±)\ / r{±)\ . /r(±) /-(±)\ r(±) 



c 



where J(±) = r c/xfsechx)^/^±+2" = ,/? ^(1/(2^^) +n) 

flx^secnxj v'^p(^/(2^2 ) + ^ + 1/2) " 

To obtain the explicit form we adopt /o = 1/(1 + tanh(2)) in (|]). In this case the 



coefficients in ([2^ ) are reduced to P23 = —3/2, P32 = — /3, p4i = —1/4, 77 = 1/(4/3) cq = 
26/0/3^ = 1.20689, and p = 3^3/(8^2/0) = 0.902037. Thus, we obtain c as 

c+ = 0.62485945; c_ = 0.82170040, (25) 
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where c± are respectively the solution of (p^ corresponding to 9±. Since there are two 
propagating velocities, the linearly stable region invades the unstable region if there are 
many domains in the system. 

To check the validity of our analysis we perform the numerical simulation of (|1]) 
and (|^) with /o = 1/(1 + tanh(2)) under the periodic boundary condition. We adopt 
the classical fourth-order Runge-Kutta scheme. Since our purpose is the quantitative 
test of (^OD and (|25|) , the initial condition is restricted to the localized symmetric form 
Tn = 18.7/A^(tanh(n — — tanh(n — 3N/A) — 1) where N is the number of cars. Tak- 
ing into account the scaling properties we perform the simulation for the set of parameters 
(e,A^) = (1/2, 32), (1/4, 64), (1/8, 128), (1/16, 256) until r„ relaxes to steady propagating 
states. Our result is plotted in Figs.l and 2. Figure 1 displays points which have the 
maximum and the minimum values of successive car distance in each parameter set, and 
theoretical coexistence curve 



« = «c(l-^^4^)' ^=\/pfj^ = 1-15850495992, (26) 

where the agreement with each other is obvious. Notice we adopt c = (c+ + c_)/2 as the 
traveling velocity, because the domain cannot move due to the finite size effect. From this 
figure we can see that one of the branches is in the linearly unstable region but the theoretical 
curve recovers the simulation result. Figure 2 demonstrates that the numerical result has 
a scaling solution which has an asymmetric kink-antikink pair. The linear combination of 
our theoretical curve (|20| ) and (p5D is plotted as the solid line by choosing the position of 
the kink and the antikink. Our theoretical curve agrees with simulation value without other 
fitting parameters. The quantitative discussion on the spreading process due to the finite 
c+ — c_ will be discussed elsewhere. 

Let us comment on the universality class of traffic flows and granular pipe flows. All of 
models introduced here except for (^) have asymmetric kink-antikink pairs and qualitatively 
resemble behaviors with each other. Komatsu |T^ has derived (|T^) as the long wave equation 



from the fluid model of traffic flow |^ which is equivalent to the two-fluid models in granular 
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flows. It is also easy to derive ([T^) from (^). In this sense, granular flows and traffic flows 
compose a universality class and our discussion here essentially can be used in any models 
for traffic flows and granular pipe flows. On the other hand, OV model in is a special case 
of the above generalized models. For example, the models with /o = in and g"{ac) = 
in (^) are reduced to MKdV equation. 

In conclusion, we obtain the analytic scaling solution of (|l]) and (|^) which describes the 
phase separation between a hnearly unstable phase and a stable phase. The accuracy and 
relevancy of the solution are confirmed by the direct simulation. 

One of the authors (HH) thanks Ooshida, T. for fruitful discussion. This work is par- 
tially supported by Grant-in-Aid of Ministry of Education, Science and Culture of Japan 
(09740314). 
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FIGURES 

FIG. 1. Plots of the coexistence curve (26), the solid line, and the neutral curve (9), the dashed 
line, as the functions of h. The scattered points arc the maximum and the minimum distances of 
successive cars for a given a obtained from our simulation. 

FIG. 2. The linear combination of our theoretical curve (20) and the scaled simulation data 
of the relative distance of succesive cars for {e,N) = (1/2, 32), (1/4, 64), (1/8, 128), (1/16, 256), 
where 'N.s' denotes the scaled data for A'^ cars systems. The solid curve is 
f{z) = tanh(^6'+(2 - z+)) - 1 +tanh(^6'_(z - z^)) with i = {Qc/cq)^/'^ /IQ = 0.11851533 and 
two fitting parameters z+ = 62.5 and z- = 190.5, where the spatial scale is measured by the 
average distance in N = 256. 
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